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Abstract 

We study in next-to-leading order QCD hadronic penguin B decays in the 
Standard and two-Higgs-doublet models. Although the gluonic penguin dom- 
inates, we find the electroweak contribution non-negligible. In the Standard 
Model, the branching ratio for B — ► X s 4> is predicted to be in the range 
(0.6 ~ 2) x 1(T 4 . The ranges of branching ratios for B — > K<f>, B — ► K*<p, 
and B s -> # are (0.4 ~ 2) x 10" 5 , (0.2 ~ 1) x 10~ 5 , and (0.15 ~ 0.5) x 10~ 5 , 
respectively. The contribution from the charged Higgs boson in two Higgs 
doublet models depend on cot(3, and can be as large as 40%. 
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Rare B decays, particularly pure penguin decays, have been subject of considerable 
theoretical and experimental interest recently [|TJ. The photonic penguin induced process 
B — > A"*7 has been observed by CLEO collaboration and is consistent with the Standard 
Model (SM) prediction 0. The gluonic penguin induced B decays are expected to be 
observed very soon. A large number of gluonic penguin induced B decay channels were 
studied in Ref. Q using AB = 1 effective Hamiltonian Hab=i in the lowest nonvanishing 
order. In Ref. || the next-to-leading order QCD corrected pure gluonic penguin -£Tab=i was 
used with top quark mass m t fixed at 150 GeV. In this paper we study the next-to-leading 
order QCD corrected Hamiltonian H&b=i in the SM and in two Higgs doublet models, 
taking particular care to include the full electroweak contributions and find the dependence 
on m t and a s . Using this Hamiltonian we study the cleanest signature of hadronic penguin 
processes: B — > X s <f), B — > K(p(K*<j)), and B s — > 00. The process B — > X s is particularly 
recommended because it is free from form factor uncertainties. We find not only that the 
QCD correction in next-to- leading order are large, but also inclusion of the full electroweak 
contributions have significant effect on the branching ratio which could reduce the pure 
gluonic penguin contribution by 30% at the upper range of allowed top quark mass. Our 
results which have been derived independently, agree with Ref. |J where only the SM is 
considered. 

AB = 1 gluonic penguin Hamiltonian 

The QCD corrected H& B= i relevant to us can be written as follows 0]: 

#ab=i = %K&C(ciC>i + c 2 2 u ) + VdVZiaOt + c 2 O c 2 ) - V tb V* J2 CiOt] + H.C. , (1) 

where the Wilson coefficients (WCs) q are defined at the scale of yU ~ and Oi are defined 

as 

0\ = SaT/x(l - Ib^pq^il - l 5 )q a , 0\ = S7 M (1 - 7 5 )6g7 M (l - 7 5 )g , 

3 = S7 M (1 -75)^^(1 -75)9' , Q4 = s a 7 M (l -75)^X]^7m( 1 -J5)q' a , 

q' q> 

5 = s llx (l - 75 )& E + 75)g' , = s al ,(l - l5 )b p E ^ 7 „(1 + l5 )q' a , (2) 
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q< q' 

3 3 
°9 = o^ni 1 ~ 75)&53 e 9 /gV(l - 7sW , Q10 = 2^7m(1 - 7s)&/3 e i'$plm u ( l ~ 7sK ■ 

q> q> 

The WCs Cj are obtained by solving the renormalization group equation 

(/^ + /^)^)C(m^V) =f(g 2 )C(m 2 w /» 2 ,g 2 ) . (3) 
Here C is the column vector (cA, and 



W) = 7<°>^ + 7i"^ + 7f £ + 7<^ + ... , (4) 

where n/ is the number of active quark flavours. 

The anomalous-dimension matrix 7^ and the first term in /3(g) determine the leading 
log QCD corrections ||. The rest of the terms contain information about the leading QED 
and next-to-leading order QCD corrections. The full 10 x 10 matrices for 7, are given in 
Ref. 0. The matching conditions of the Wilson coefficients at raw for the next-to- leading 
order corrections will be different from the leading order ones. One needs to include one 
loop current-current corrections for at raw- The full results for the initial conditions can 
be found in [0. 

The WCs obtained above depend on the renormalization regularization scheme (RS) 
used. In our calculation we used the naive dimentional regularization scheme. The physical 
quantities, of course, should not depend on RS provided one handels the hadronic matrix 
elements correctly. In practice, many of the hadronic matrix elements can only be calculated 
using factorization method. In our later calculation we will also use this approxmation. Since 
this approximation does not carry information about the RS dependence, it is better for us 
to use WCs, C(/i) = (1 + fjfa s (/x)/47T + fgQ; em (^)/47r)C(//), which are RS independent ||. 
Here the matrices f s e are obtained from one-loop matching conditions. The 6 x 6 f s matrix 
for the pure gluonic penguin operators has been given in Ref. ||. Based on the work of Ref. 
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||, we have worked out the full 10 x 10 matrices for r Sje and carried out the calculation 
using the full matrices. 

We also need to treat the matrix elements to one-loop leve for consistency. These one- 
loop matrix elements can be rewritten in terms of the tree-level matrix elements < Oj > tree 
of the effective operators, and one obtains PJTI1|] 

< cA >= £ c,0u) [<% + + ^m%] < O >*~ . (5) 

We have worked out the full matrices m s,e . For the processes we are considering only C3_io 
contribute. Expressing the effective coefficients which multiply the matrix elements 
< Oj >* ree in terms of Cj, we have 

eff — r) /q e// — \ TD e ff — JD /Q e // — i 73 

C 3 — C3 — -Ts/O , C 4 — C3 + J^s , C 5 — C5 — r s / 6 , C 6 — C6 + r s , 

c e 7 ff = c 7 + P e , cj f = c 8 , 4 ff = c 9 + P e , cll f = C10 • (6) 

The leading contributions to P SiE are given by: P s = (a s /87r)c 2 (10/9 + G(m c , /i, q 2 )) and 
P e = (a em /97r)(3ci + C2)(10/9 + G(m c , fi,q 2 )). Here m c is the charm quark mass which we 
take to be 1.35 GeV. The function G(m,/i, q 2 ) is give by 

G(m,/i,g 2 ) = 4 / x(l - x)dxln ^ — . (7) 

In the numerical calculation, we will use q 2 = ml/2 which represents the average value and 
the full expressions for P s e . 

Using range of values of a s {mz) and m t we can calculate the coefficients at fi = nib. We 
use a s {m z ) as input instead of A 4 as in Ref. 0. In Table 1, we show some sample WCs 



for the central world average value of a s (mz) = 0.118 [|TTJ and for several values of m t with 
a em = 1/128. 

In the two-Higgs-doublet model, there are new contributions to q due to charged Higgs 
boson. The charged Higgs-quark couplings are given by [|12j 



L H = 2 J^ m Ui Vij \. cot P m ^ ( 1 ~ 75) - am dj . (1 + ry 5 )]djH + + H.C. , (8) 



where cot (3 = V\jv%\ v\ and v<i are the vacuum expectation values of the Higgs doublets 
Hi and H 2 , which generate masses for down and up quarks, respectively. The parameter a 



depends on the models (T^|. The main contributions are from the first term in eq.(8) and 
we will neglect the contribution from the second term. The charged Higgs contributions 
to gluonic penguin have been studied by several groups |T3. The leading QCD corrected 



Hamiltonian has been given in Ref. ||14j| . We have checked the next-to- leading initial con- 



ditions for the WCs at my/- We find that the inclusion of charged Higgs will not change 
the initial conditions for 01,2,8,10) but 03,4,5,6,7,9 are changed in the same way as those given 



in eqs. (32-38) of Ref. |J. 
B->X s 

Using H&b=i in eq.(l), we obtain the decay amplitude for B — > X s (p 

A(B -> X^) « A(b -> s<P) = -^ Vtb V*e»Cs lti (l - l5 )b , (9) 

where e M is the polarization of the </> particle; C = cjj + C4 + C5 + £(cj} + C4 + Cg ) — 
(cy + Cg + cj/ + £ (c§ + Cg + c\q ))/2 with £ = l/iV c , where N c is the number of colors. 
The coupling constant is defined by < (f)\s^s\Q >= ig^. From the experimental value 
for Br((j) — > e + e~), we obtain g^ = 0.0586 GeV A . 
The decay rate is, then, given by 

T(B - X s 0) = ^^^^[l + _ + 4)] , (10) 

where A^- = (1 — m 2 /m1 — m 2 /m 2 ) 2 — 4m 2 m 2 /mf. 

We normalize the branching ratio to the semi-leptonic decay of B — > X c ei> e . We have 



Sr(S -> X s 0) = 5r(B -> X c et/ e ) 5—^ C 1 + 1 (1 1 H f ) . (11) 



In the above expression, p = 1 — 8r 2 + 8r 6 — r 8 — 24r 4 lnr with r = m c / nib, is the phase factor, 
and f) is the QCD correction factor in B — > X c ez/ e , respectively. We will use p = 0.5, 77 = 
0.889 and the approximation | VtbV t */V c b\ 2 = 1. The branching ratio Br(B — > X c e^ e ) is 



measured to be 0.108 116 



Exclusive decays B — > K(j){K*<p) and £> s — > 00 

For the exclusive decays, we will use the factorization method. We have 

Gp . 



A(B -> K{K*)4>) 
A(B S - 



V2 
Gp 

"71 



< X(^*)|s 7 ^(l - lb )b\B >< 0|s 7 ^|O > 
V tb V*C < 0|^(1 - 7 5 )&|5 S >< <p\^s\0 > . 



(12) 



We can parametrize the matrix elements as 



< K\s lfl (l - l5 )b\B > = f + (q 2 )(p*+p«) + f-(q 2 )q, 



<v\s lil {\- 1 ^)b\B> = 2V(q 2 )ie, ^->*-»" 



(13) 



+ Ax(g 2 )K - ™bK ~ M{q')e v ■ q(p B + p 



where v is a vector meson particle and its polarization. For B — > Kcj), q = p — p , and 



for B — > v(p, q = p B — p v . 



In terms of the form factors defined above, we obtain the decay rates 
G 2 F f +2 (m 2 )g 2 mj 



T{B - 
T(B 



K<j>) 



IB \v th vz\ 2 c 2 \T. 



K<t> 5 



v<P) = ^^\V tb V;fC 2 X% 2 [2V 2 (m 2 ) + ±{1 - ^IH) - A\(m\ 
1 



+ 



V(f> 

2 _2 



Am 2 



v m '<i> 



-(« - m^)Ai(mJ) - Kj - m 2 - m 2 )A 2 (m 2 )) 2 ] 



(14) 



To finally obtain the branching ratios, we will use two sets of form factors obtained by 
Bauer et. al. |T?| and Casalbuoni et. al. fll~8| . Note that we have used different normalization 



for the form factors V and A4 from those in Refs. 1 17,18]. The form factors at q = are 



determined by using relativistic quark model in Ref. JT^, and by using chiral and effective 
heavy quark theory in Ref. |l£|. The form factors at q 2 = in Ref. [17| are given by: 
/+^ = 0.393, V K *^ = 0.062 GeV~\ A lK «^ = -0.077 GeV~\ and A 2K , 4> = 0.056 GeV' 1 . 
In Ref. the form factors at q 2 = are: = 0.509, V K *$ = 0.103 GeV~\ A^*^ = 
-0.047 GeV~\ and A 2K *<t> = 0.034 GeV' 1 . In Ref. the form factors at g 2 = for 
B s — > 00 are also calculated. They are: = 0.102 GeV' 1 , A\^ = —0.046 GeV -1 , and 
Ai^ = 0.033 GeV" 1 . In both papers, the q 2 dependence of all the form factors were assumed 
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to be of a simple pole type. We will use the pole masses used in Refs. [|T^,|T^]. It is interesting 
to note that the ratios between the exclusive decays and B — > X s <ft are independent of the 
Wilson coefficients. If these ratios can be measured experimentally, they can test the models 
for the form factors. We obtain 

0.06 , Ref. [17] 
0.1 , Ref. g 



(15) 



Br(B - 




Br(B - 


+ x s <p) 


Br(B - 


>K*<j>) 


Br(B - 


>X s ci>) 


Br(B - 


-00) 



0.057 , Ref. 



0.029 , Ref. [H3 



0.023 , Ref. HI 



Br(B -> X s <Pj 

We show in, Fuigure 1, the predictions for the branching ratio Br(B — > X s <p) in the 
SM as a function of top quark mass m t and the strong coupling constant a s (mz). The 
QCD corrections turn out to be important which enhance the branching ratios by about 
30% compared with those of without QCD corrections. There is a large uncertainty in the 
branching ratios due to error in a s (mz). From Figure 1, we see that the error in a s (mz) 
can induce an uncertainty of a factor 2. 

The dominant contribuitons are from the gluonic penguin. There is a very small m t 
dependence for the branching ratio calculated without the inclusion of the electroweak pen- 
guin contributions. The inclusion of the full electroweak contribuitons have sizeable effects 
which reduce the branching ratios by about 20% to 30% for the central value of a s with m t 
varying from 100 GeV to 200 GeV. It is clear from Figure 1 that the full contribution has a 
large m t dependence. 

There may be corrections to the branching ratios predicted by the factorization method. 
It is a common practice to parameterize the possible new contributions by treating £ as a 
free parameter ||i7|-|T9"| . Using experimental values from non-leptonic B decays, it is found 



that ||18||, a\ = C2 + £ci and 02 = Ci+£c2 have the same signs, and [0.2 1 ~ 0.27 and |ai| ~ 1.0. 



We see that £ is close to 1/2. To see the effect of varying £, we plot the predictions for the 



branching ratios for £ = 1/2 and £ = 1/3. The branching ratios for £ = 1/2 are about 2 
times those for £ = 1/3. 

For the central value of a s (mz) and the central value of m t = 174 GeV reported by 
CDF |2(J, the value for Br(B -> V s 0) is about 1.7 x 1(T 4 for £ = 1/2. The exclusive 
branching ratios 5 — > X0 and I? — > ii"*0 are about the same which are 1 x 10~ 5 if the 
form factors from Ref. Jl7| are used. If the form factors from Ref. |l8j are used, one obtains 
Br(B -> K(j>) « 1.7 x 10" 5 , £r(£ -> K*0) « 0.5 x 10~ 5 , and Br(B s -> 00) « 0.4 x 10~ 5 . 

In Figure 2, we show the ratio of the branching ratios BriB — > X s <f))2H and Br(B — > 
X s <P)sm predicted by the two Higgs doublet model and the SM as a function of cot/3 for 
mt = 174 GeV and different values of m# with £ = 1/2. The depence on £ is small. From 
Figure 2, we see that the effects of the charged Higgs boson contributions are small for 
cot fi < 1. When increasing cotfi, the charged Higgs contributions become important and 
the effect is to cancel the SM contributions. When cotfi becomes very large the charged 
Higgs boson contributions become the dominant ones. However, using the information from 
B — > X s 7, it is found that for small m# ~ 100 GeV and m t ~ 174 GeV, cotfi is constrained 
to be less than 1 |f2TJ. For these values, the charged Higgs boson effects on the processes 
discussed in this paper are less than 10%. For m H ~ 500 GeV, the charged Higgs boson 
effects can reduce the hadronic penguin B decays by 40% because the range of cotfi allowed 



from b — > 57 is now larger |[21|| . The effects become smaller for larger m#. 

The analyses carried out in this letter can be generalized to other hadronic B decays. 
We will present the full calculations for the Wilson coefflceints, the full expressions for P s e 
and other related decays in a forthcoming paper |22| . 



We thank Buras, McKellar, Fleischer for useful corespondences and thank Lautenbacher 
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TABLES 



TABLE I. The Wilson coefficients for AB = 1 at m b = 5 GeV with a s {m z ) = 0.118. 



mt(GeV) 


Cl 


C2 


C3 


c 4 


C5 


Q3 


C7/a em 


C8/a em 




Cio/Oem 


130 


-0.313 


1.150 


0.017 


-0.037 


0.010 


-0.045 


-0.061 


0.029 


-0.978 


0.191 


174 


-0.313 


1.150 


0.017 


-0.037 


0.010 


-0.046 


-0.001 


0.049 


-1.321 


0.267 


210 


-0.312 


1.150 


0.018 


-0.038 


0.010 


-0.046 


0.060 


0.069 


-1.626 


0.334 
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Figure Captions 



Figure 1. Br(B — > X s (p) as a function of m t and a s (m z ). The regions between the 
dashed and solid lines are the branching ratios for a s {mz) varying from 0.111 to 0.125 for 
£ = 1/2 and £ = 1/3, respectively. The branching ratios increases with a s (mz)- 

Figure 2. Br(B — > X s 0) 2 h / 1 Br(B — > X s 4>)sm as a function of cot/3, and m^. The curves 
1, 2 and 3 are for m H equals to 100, 500 and 1000 GeV, respectively. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9403266v3 



